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ABSTRACT 

We  have  assembled  a  spectroscopic  sample  of  low-mass  dwarfs  observed  as  part  of  the  Sloan  Digital  Sky  Survey 
along  one  Galactic  sight  line,  designed  to  investigate  the  observable  properties  of  the  thin  and  thick  disks.  This  sample 
of  ^7400  K  and  M  stars  also  has  measured  ugriz  photometry,  proper  motions,  and  radial  velocities.  We  compute  UVW 
space-motion  distributions,  and  investigate  their  structure  with  respect  to  vertical  distance  from  the  Galactic  plane. 
We  place  constraints  on  the  velocity  dispersions  of  the  thin  and  thick  disks,  using  two-component  Gaussian  fits.  We 
also  compare  these  kinematic  distributions  to  a  leading  Galactic  model.  Finally,  we  investigate  other  possible  observ¬ 
able  differences  between  the  thin  and  thick  disks,  such  as  color,  active  fraction,  and  metallicity. 

Key  words:  Galaxy:  kinematics  and  dynamics  —  Galaxy:  structure  —  stars:  fundamental  parameters  — 
stars:  low-mass,  brown  dwarfs 


1.  INTRODUCTION 

Modeling  the  Galaxy  presents  a  challenging  breadth  of  prob¬ 
lems  to  both  theorists  and  observers.  Large  X-body  simulations 
employing  the  constraints  of  gravity  and  A  cold  dark  matter  cos¬ 
mology  have  sought  to  recreate  the  infant  Galaxy,  tracing  the 
formation  and  collapse  of  baryons  within  the  dark  matter  halo 
(Brook  et  al.  2004;  Governato  et  al.  2007).  Observationally, 
these  simulations  are  constrained  by  rotational  velocities  and 
luminosity  profiles  of  extragalactic  systems  ( Dalcanton  et  al. 
1997).  Closer  to  home,  observers  seek  to  reconstruct  the  his¬ 
tory  of  the  Milky  Way’s  cannibalistic  mergers  through  photo¬ 
metric  identification  of  tidal  debris,  such  as  the  Sagittarius  dwarf 
(Ibata  et  al.  1994).  Spectroscopy  is  also  employed  to  find  sur¬ 
viving,  comoving  stars  of  similar  metallicity  (e.g.,  Yanny  et  al. 
2003). 

A  particularly  interesting  problem  for  both  the  theorist  and 
observer  is  the  formation  and  nature  of  the  thick  disk  (Gilmore 
&  Reid  1983).  This  population  has  been  explored  extensively, 
mainly  through  star  counts  (e.g.,  Gilmore  &  Reid  1983;  Reid  & 
Majewski  1993;  Buser  et  al.  1999;  Norris  1999;  Siegel  et  al. 
2002).  However,  the  thick  disk  scale  height  and  local  density 
normalization  are  still  uncertain  (Norris  1999),  with  values  of 
^ thick  ranging  from  ^700  to  1500  pc  and  local  normalizations 
between  2%  and  15%.  Larger  scale  heights  are  usually  coupled 
to  lower  normalization  values  (see  Fig.  1  of  Siegel  et  al.  2002). 
The  scale  length  of  the  thick  disk  is  also  uncertain,  although 
typically  scale  lengths  larger  than  the  thin  disk  are  inferred  (Chen 
et  al.  2001;  Larsen  &  Humphreys  2003).  The  thick  disk  is 
thought  to  be  an  older,  metal-poor  population  ( Reid  &  Majewski 
1993;  Chiba  &  Beers  2000).  Metallicity  differences,  such  as 
a -element  enhancement,  have  been  explored  (Bensby  et  al. 
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2003),  but  other  observables,  such  as  photometric  color  and 
chromospheric  activity  trends,  have  yet  to  be  studied  in  detail. 
Thick  disks  observed  in  external  galaxies  ( Yoachim  &  Dalcanton 
2006)  often  appear  to  have  structural  parameters  and  kinematics 
similar  to  the  Milky  Way. 

Modeling  the  observable  properties  of  the  thin  and  thick  disks, 
namely,  star  counts  (Reid  &  Majewski  1993;  Reid  1993)  and 
kinematics  (Mendez  &  van  Altena  1996;  Robin  et  al.  2003; 
Vallenari  et  al.  2006),  has  undergone  a  resurgence  in  recent 
years.  The  seminal  work  of  Bahcall  &  Soneira  (1980)  laid  the 
foundation  for  modeling  star  counts  of  the  smooth  Galactic 
components.  Present-day  models,  such  as  the  Besamjon  model 
(Robin  et  al.  2003)  and  the  Padova  model  (Vallenari  et  al.  2006), 
also  incorporate  kinematics,  allowing  for  robust  comparisons  to 
observations.  It  is  important  to  rigorously  test  their  predictions 
against  actual  observed  spatial  and  kinematic  distributions. 

Low-mass  dwarfs  are  both  ubiquitous  and  long-lived  (Laughlin 
et  al.  1997),  and  serve  as  excellent  tracers  of  the  Galactic  po¬ 
tential  (Wilson  &  Woolley  1970;  Wielen  1977;  Westetal.  2006). 
Modem  surveys,  such  as  the  Sloan  Digital  Sky  Survey  (SDSS; 
York  et  al.  2000),  are  sensitive  to  K  and  early  M  dwarfs  at  dis¬ 
tances  of  ~  1  kpc  above  the  Galactic  plane,  probing  the  transition 
between  the  thin  (Reid  et  al.  1997)  and  thick  (Kerber  et  al.  2001) 
disks.  At  these  distances,  we  expect  about  20%  of  the  observed 
stars  to  be  thick  disk  members,  assuming  a  local  normalization 
of  2%  (Reid  &  Majewski  1993)  and  scale  heights  of  300  and 
1400  pc  for  the  thin  and  thick  disks,  respectively.  Star  counts  of 
low-mass  dwarfs  have  been  used  to  determine  Galactic  structural 
properties  (Reid  et  al.  1997  and  references  therein).  These  stud¬ 
ies  sought  to  determine  the  vertical  scale  height  and  local  nor¬ 
malization  of  the  thin  and  thick  disks  (Siegel  et  al.  2002),  as  well 
as  the  underlying  mass  function  of  these  populations  (Martini 
&  Osmer  1998;  Phleps  et  al.  2000;  K.  C.  Covey  et  al.  2008,  in 
preparation).  The  chemical  evolution  of  the  Galaxy  has  also 
been  explored  with  red  dwarfs  (Reid  et  al.  1997),  using  mo¬ 
lecular  band  indices  as  a  proxy  for  metallicity  (Gizis  1997). 

In  addition  to  their  utility  as  Galactic  tracers,  low-mass  dwarfs 
have  intrinsic  properties,  such  as  chromospheric  activity  (West 
et  al.  2004;  Bochanski  et  al.  2005;  Schmidt  et  al.  2007),  that  can 
be  placed  in  a  Galactic  context.  In  addition,  metallicity  may  be 
probed  using  subdwarfs  (Lepine  et  al.  2003),  readily  identified 
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TABLE  1 

Previous  Low-Mass  Stellar  Kinematic  Survey  Results 


Study 

Sample  Size 

Distance  Limit 
(pc) 

a 

GV 

aw 

Comments 

Vyssotsky  (1956);  Reid  et  al.  (1995) . 

368 

~300 

39 

26 

23 

McCormick  stars  in  PMSU  survey 

Wielen  (1977) . 

516 

~300 

39 

23 

20 

McCormick  stars 

Reid  et  al.  (1995) . 

514 

25 

43 

31 

25 

PMSU  I  volume-complete  sample 

Ratnatunga  &  Upgren  (1997) . 

773 

~300 

30.6 

18.5 

7.4 

McCormick  stars 

Reid  et  al.  (2002) . 

436 

25 

37.9 

26.1 

20.5 

PMSU  IV  volume-complete  sample 

25 

34 

18 

16 

PMSU  IV  volume-complete  sample,  core 

Bochanski  et  al.  (2005) . 

419 

100 

35 

21 

22 

Nonactive  stars,  core 

155 

100 

19 

15 

16 

Active  stars,  core 

This  study . 

70 

100 

28.4 

21.3 

19.2 

|z|  <  100  pc 

100 

25.7 

20.9 

14.1 

|z|  <100  pc,  core 

4300 

500 

39.0 

30.0 

24.8 

\z\  <  500  pc 

500 

34.1 

24.6 

20.9 

|z|  <  500  pc,  core 

6893 

1000 

44.1 

36.9 

27.9 

|z|  <  1000  pc 

1000 

37.4 

27.2 

23.6 

|z|  <  1000  pc,  core 

7398 

1600 

46.6 

42.2 

29.3 

Total  sample 

1600 

38.5 

28.6 

24.4 

Total  sample,  core 

The  velocity  dispersions  in  U,  V,  and  W  are  measured  in  km  s  1 


by  their  spectra,  which  show  enhanced  calcium  hydride  (CaH) 
absorption.  Subdwarfs  have  been  easily  detected  in  large-scale 
surveys  such  as  the  SDSS  (West  et  al.  2004). 

Kinematic  studies  of  low-mass  stars  have  a  rich  historical  back¬ 
ground.  Samples  are  typically  drawn  from  proper-motion  surveys, 
such  as  the  New  Luyten  Two  Tenths  (NLTT)  catalog  (Luyten 
1979)  and  the  Lowell  Proper  Motion  Survey  (Giclas  etal.  1971). 
Efforts  have  been  made  to  identify  nearby  stars  in  these  surveys, 
for  example,  by  Gliese  &  Jahreiss  (1991 ).  However,  proper-motion 
surveys  possess  inherent  kinematic  bias.  The  McCormick  sam¬ 
ple  (Vyssotsky  1956),  assembled  from  875  spectroscopically 
confirmed  K  and  M  dwarfs,  has  been  frequently  studied  as  a 
kinematically  unbiased  sample  (Wielen  1977;  Weis  &  Upgren 
1995;  Ratnatunga  &  Upgren  1997),  although  it  has  been  sug¬ 
gested  that  the  sample  may  be  biased  toward  higher  space  mo¬ 
tions  (Reid  et  al.  1995).  The  Palomar-Michigan  State  University 
survey  (PMSU;  Reid  et  al.  1995, 2002;  Hawley  et  al.  1996;  Gizis 
et  al.  2002)  is  among  the  largest  prior  spectroscopic  surveys  of 
low-mass  stars,  obtaining  spectral  types  and  radial  velocities  of 
~1700  M  dwarfs.  The  PMSU  sample,  which  targeted  objects 
from  the  Third  Catalogue  of  Nearby  Stars  (Gliese  &  Jahreiss 
1991  ),  was  used  to  construct  a  volume-complete,  kinematically 
unbiased  catalog  of  ~500  stars,  sampling  distances  to  ~25  pc. 
Other  surveys,  such  as  the  100  pc  survey  (Bochanski  etal.  2005), 
have  also  used  low-mass  stars  as  kinematic  probes.  In  Table  1  we 
summarize  the  sample  sizes  and  approximate  distance  limits  of 
previous  major  kinematic  surveys  of  low-mass  dwarfs,  along 
with  the  mean  velocity  dispersions  for  each  study.  It  is  clear  that 
our  sample  of  several  thousand  stars  out  to  distances  of  ~1  kpc 
results  in  a  study  of  Galactic  kinematics  using  low-mass  stars 
with  unprecedented  statistical  significance. 

In  this  paper  we  present  our  examination  of  the  properties 
of  the  thin  and  thick  disks  using  an  SDSS  Low-Mass  Spectro¬ 
scopic  Sample  (SLoMaSS)  of  K  and  M  dwarfs  that  is  an  order  of 
magnitude  larger  than  previous  samples.  In  §  2  we  describe  the 
SDSS  spectroscopic  and  photometric  observations  that  comprise 
SLoMaSS.  The  resulting  distances  and  stellar  velocities  are  pre¬ 
sented  in  §  3.  In  §  4  we  detail  our  efforts  to  separate  the  observed 
stars  into  two  populations,  search  for  kinematic,  metallicity,  and 
color  gradients,  and  compare  our  results  to  a  contemporary  Ga¬ 
lactic  model.  Finally,  §  5  summarizes  our  findings. 


2.  OBSERVATIONS 

2.1.  SDSS  Photometry 

The  SDSS  (York  et  al.  2000;  Stoughton  et  al.  2002;  Pier  et  al. 
2003;  Ivezic  et  al.  2004)  is  a  large  (~  10,000  deg2),  multicolor 
(ugriz;  Fukugita  et  al.  1996;  Gunn  etal.  1998;  Hogg  et  al.  2001; 
Smith  et  al.  2002;  Tucker  et  al.  2006)  photometric  and  spec¬ 
troscopic  survey  centered  on  the  northern  Galactic  cap.  The 
2.5  m  telescope  (Gunn  et  al.  2006),  located  at  Apache  Point  Ob¬ 
servatory,  scans  the  sky  on  great  circles,  as  the  camera  (Gunn  et  al. 
1998)  simultaneously  images  the  sky  in  five  bands  to  a  faint 
limit  of  ~22.2  in  r,  with  a  typical  uncertainty  of  ~ 2 %  at  r  ~  20 
(Ivezic  et  al.  2003;  Adelman-McCarthy  et  al.  2006).  The  last 
data  release  (DR5;6  Adelman-McCarthy  et  al.  2007)  comprises 
8000  deg2  of  imaging,  yielding  photometry  of  ~217  million 
unique  objects,  including  ~85  million  stars.  SDSS  photom¬ 
etry  has  enabled  a  myriad  of  studies  on  both  Galactic  structure 
(e.g.,  Yanny  et  al.  2000;  Newberg  et  al.  2002;  Juric  et  al.  2005; 
Belokurov  et  al.  2006)  and  low-mass  stars  (e.g.,  Hawley  et  al. 
2002;  Walkowicz  et  al.  2004;  West  et  al.  2005;  Davenport  etal. 
2006). 

2.2.  SDSS  Spectroscopy 

When  sky  conditions  prohibit  photometric  observations, 
the  SDSS  telescope  is  fitted  with  twin  fiber-fed  spectrographs. 
These  instruments  simultaneously  obtain  640  medium-resolution 
(R  ~  1 800),  flux-calibrated,  optical  (3800-9200  A)  spectra  per  3° 
plate,  permitting  radial  velocity  measurements  for  most  stars  with 
an  uncertainty  of  ~  10  km  s~*  (Abazajian  et  al.  2004).  A  typi¬ 
cal  45  minute  observation  yields  a  signal-to-noise  ratio  per  pixel 
greater  than  4  at  g  =  20.2  and  i  =  19.9  (Stoughton  et  al.  2002), 
with  a  broadband  flux  calibration  uncertainty  of  ~ 4 %  (Abazajian 
et  al.  2004).  The  DR5  sample  includes  over  1  million  spectra, 
with  ~2 16,000  stellar  spectra  (  Adelman-McCarthy  et  al.  2007). 
The  majority  of  spectra  in  the  SDSS  database  are  drawn  from  three 
main  samples,  which  target  objects  based  on  their  photometric 
colors  and  morphological  properties.  These  samples  are  opti¬ 
mized  to  observe  galaxies  (Strauss  et  al.  2002),  luminous  red  gal¬ 
axies  with  r  ~  0. 5-1.0  (Eisenstein  et  al.  2001),  and  high-redshift 


6  See  http://www.sdss.org/dr5/. 
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Fig.  1. — Top :  Distribution  of  spectral  types  within  SLoMaSS.  The  mean  i  —  z  color  at  each  spectral  type  is  shown  above  each  bin.  Bottom :  The  i  vs.  i  —  z  color- 
magnitude  distribution  and  r  —  i  vs.  i  —  z  color-color  diagram  of  stars  in  SLoMaSS.  The  lowest  contour  (black)  in  the  color-color  diagram  indicates  a  density  of  100  stars 
in  a  square  0. 1  mag  wide.  Each  additional  contour  represents  100  additional  stars. 


quasars  ( Richards  et  al.  2002).  However,  low-mass  stars  have  sim¬ 
ilar  colors  to  some  of  these  samples  and  therefore  are  observed 
serendipitously.  SDSS  spectroscopy  of  late-type  dwarfs  has  been 
the  focus  of  numerous  studies  (see  Hawley  et  al.  2002;  Raymond 
et  al.  2003;  West  et  al.  2004, 2006;  Silvestri  etal.  2006;  Bochanski 
et  al.  2007  and  references  therein). 

2.3.  SDSS  Low-Mass  Spectroscopic  Sample:  SLoMaSS 

During  fall  2001,  a  call  was  placed  to  the  SDSS  collaboration 
to  design  special  spectroscopic  plates  that  employed  different 
targeting  algorithms  than  the  usual  SDSS  survey  samples  de¬ 
scribed  in  §  2.2.  We  designed  and  proposed  a  series  of  obser¬ 
vations  to  probe  the  local  vertical  structure  of  the  Milky  Way, 
obtaining  spectra  of  low-mass  dwarfs  in  the  southern  equatorial 
stripe  (stripe  82)  of  the  SDSS  photometric  footprint  during  fall 
2002.  This  stripe  at  zero  declination  is  repeatedly  observed  dur¬ 
ing  the  time  of  the  year  when  the  northern  Galactic  cap  is  not 
visible.  These  repeat  scans  sample  an  area  of  ^300  deg2  and 
have  been  used  to  study  Type  la  supernovae  (Sako  et  al.  2005; 
Frieman  et  al.  2007)  and  stellar  variability  (Sesar  et  al.  2007) 
and  characterize  the  repeatability  of  SDSS  photometry  (Tvezic 
et  al.  2007).  SLoMaSS  comprises  stars  with  unsaturated  griz 
photometry  and  extinction-corrected  magnitude  limits  of  1 5  < 
i  <  18  and  i  —  z  >  0.2.  A  series  of  three  spectroscopic  tilings 
composed  of  15  plates  (numbers  1118-1132,  centered  on  /  ~ 
105°,  b  ~  —62°)  was  observed,  yielding  a  total  of  8880  stellar 
spectroscopic  targets  in  SLoMaSS,  each  with  ugriz  photometry. 


3.  ANALYSIS 
3.1.  Spectral  Types 

The  spectroscopic  data  were  analyzed  with  the  HAMMER 
suite  of  software  (Covey  et  al.  2007),  which  measures  spectral 
type,  Ha  emission-line  properties,  and  various  spectral  band 
indices.  This  pipeline  automatically  assigns  a  spectral  type  to 
each  star,  then  allows  the  user  to  confirm  and  adjust  the  spectral 
type  manually,  if  necessary.  The  accuracy  of  the  final  spectral 
types  is  ±  1  subtype.  After  examining  the  entire  sample  by  eye, 
we  selected  spectroscopically  confirmed  K  and  M  dwarfs.  This 
cut  decreased  the  sample  from  8880  to  8696  stars.  The  final  spec¬ 
tral  type  distribution  of  SLoMaSS,  after  the  additional  cuts  ex¬ 
plained  below,  together  with  the  average  i  —  z  color  for  each 
spectral  type  bin,  is  shown  in  Figure  1. 

3.2.  Distances:  Photometric  Parallax 

Distances  to  the  SLoMaSS  stars  were  determined  using  the 
photometric  parallax  relations  described  in  West  et  al.  (2005)  and 
Davenport  et  al.  (2006).  The  West  et  al.  (2005)  relation  was  em¬ 
ployed  for  stars  with  i  —  z  redder  than  0.37,  and  the  Davenport 
et  al.  (2006)  parallax  relation  was  applied  to  the  bluer  stars  in 
SLoMaSS.  We  neglect  reddening  corrections,  since  the  average 
extinction  computed  for  the  total  column  along  this  line  of  sight 
from  Schlegel  et  al.  (1998)  was  small  (<0.05  in  i),  and  most  of 
these  stars  will  lie  in  front  of  significant  amounts  of  dust.  Stars 
with  i  —  z  colors  that  did  not  fall  within  the  appropriate  boundaries 
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Fig.  2. — Absolute  vertical  distance  distribution  of  the  sample  binned  every 
100  pc.  The  solid  thick  histogram  is  the  total  distribution,  while  the  five  remain¬ 
ing  histograms  represent  the  height  distribution  for  five  i  —  z  color  bins  as  de¬ 
scribed  in  the  legend.  Note  that  there  are  hundreds  of  stars  per  bin  out  to  an 
absolute  vertical  distance  of  ~1000  pc. 

of  the  photometric  parallax  relations  (0.22  <  i  —  z  <  1.84)  were 
removed  from  the  sample,  decreasing  its  size  to  8280  stars.  In 
addition,  we  applied  the  photometric  white  dwarf — M  dwarf 
pair  cuts  of  Smolcic  et  al.  (2004),  removing  four  more  stars 
from  the  sample.  The  i,  i  —  :  Hess  diagram  and  r  —  i,i  —  z  color- 
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color  diagram  are  shown  in  Figure  1  (bottom).  The  distance  from 
the  Galactic  plane  was  computed  assuming  the  Sun’s  vertical 
position  to  be  15  pc  above  the  plane  (Cohen  1995;  Ng  et  al. 
1997;  Binney  et  al.  1997).  The  vertical  distribution  of  stars  in 
SLoMaSS  is  shown  in  Figure  2.  Note  that  the  decline  in  stellar 
density  at  large  distances  reflects  the  incompleteness  of  our  sam¬ 
ple,  not  the  intrinsic  stellar  density,  while  the  decline  at  small  dis¬ 
tances  from  the  plane  is  due  to  the  saturation  of  SDSS  photometry. 

3.3.  Distances-Spectroscopic  Parallax 

In  addition  to  the  photometric  parallax  relations  mentioned 
above,  we  also  computed  the  distances  to  the  M  dwarfs  in 
SLoMaSS  using  the  spectroscopic  parallax  relations  of  Hawley 
et  al.  (2002).  These  distances  served  as  a  check  on  the  photo¬ 
metric  parallax  relations  and  were  later  used  to  divide  the  sample 
in  a  search  for  color  differences  between  the  thin  and  thick  disks 
(see  §  4.3.2  for  details).  The  spectroscopic  parallax  of  K  stars  was 
not  computed,  as  no  reliable  calibrated  relations  were  available. 

3.4.  Radial  Velocities,  Proper  Motions,  and  UVW  Velocities 

Radial  velocities  were  computed  by  cross-correlating  the 
M  dwarf  stellar  spectra  against  the  low-mass  template  spec¬ 
tra  of  Bochanski  et  al.  (2007)  with  the  IRAF7  task  fxcor.  This 


7  IRAF  is  distributed  by  the  National  Optical  Astronomy  Observatory,  which  is 
operated  by  the  Association  of  Universities  for  Research  in  Astronomy,  Inc.,  under 
cooperative  agreement  with  the  National  Science  Foundation. 
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Fig.  3. — Schematics  of  velocity  distributions  (left)  and  their  corresponding  probability  plots  (right).  The  top  row  displays  two  Gaussian  distributions,  one  with  a 
standard  deviation  of  ~50  km  s-1  (solid  line)  and  another  with  a  standard  deviation  of  ~  100  km  s-1  (dashed  line).  A  straight  line  in  the  probability  plot  signifies  an 
underlying  Gaussian  distribution,  and  the  slope  of  the  line  is  a  measure  of  the  standard  deviation  of  the  distribution.  Note  that  the  slope  of  the  dashed  line  is  twice  that  of 
the  solid  line.  In  the  bottom  row,  a  low-dispersion  component  (dashed  line)  and  a  high-dispersion  component  (dotted  line)  are  summed  to  produce  the  solid  histogram. 
While  the  dual  nature  of  the  distribution  is  not  readily  apparent  in  the  solid  histogram,  it  is  easily  detected  in  the  probability  plot. 
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Fig.  4. — Illustrative  example  of  our  velocity  dispersion  analysis.  Shown  are  probability  plots  for  U,  V,  and  W  for  the  1200  pc  <  z  <  1300  pc  bin  (solid  lines).  The  core 
component  (dashed  lines)  is  fit  between  ±  1  a,  while  the  wing  component  (dot-dashed  lines)  fits  the  outer  edges  of  the  distribution.  Note  the  strong  asymmetric  drift 
component  manifested  as  a  change  in  slope  at  negative  V  velocities. 


cross-correlation  routine  is  based  on  the  method  described  in 
Tonry  &  Davis  (1979).  The  measured  velocities  have  external 
errors  of  ~  4  km  s_1.  The  K  dwarf  radial  velocities  were  ob¬ 
tained  directly  from  the  SDSS/Princeton  ID  spectral  pipeline,8 
with  typical  uncertainties  of  10  km  s_1. 

Proper  motions  were  determined  from  the  SDSS+USNO-B 
proper-motion  catalog  (Munn  et  al.  2004).  The  catalog  is  ~90% 
complete  over  the  magnitude  limits  of  our  sample,  with  random 
errors  of  ~3. 5  masyr-1.  We  removed  an  additional  878  stars  with 
poorly  measured  proper  motions,9  reducing  our  final  sample  size 
to  7398  stars. 

Combining  distances,  proper  motions,  and  radial  velocities, 
we  computed  the  UVW  space  motions  of  each  star  in  SLoMaSS 
using  the  method  of  Johnson  &  Soderblom  (1987).  The  veloc¬ 
ities  are  computed  in  a  right-handed  coordinate  system,  with 
positive  U  velocity  directed  toward  the  Galactic  center  and 
corrected  for  the  solar  motion  (10,  5,  and  7  km  s  ,  Dehnen  & 
Binney,  1998)  with  respect  to  the  local  standard  of  rest. 

8  See  http://spectro.princeton.edu. 

9  Specifically,  we  required  a  match  between  catalogs  (match  >  0),  detections 
in  at  least  four  plates  in  USNO-B  (nfit  >  5),  no  other  objects  within  7",  which  is 
the  resolution  of  the  POSS  plates  (dist22  >  7;  see  Kilic  et  al.  2006),  and  small 
errors  in  the  proper-motion  determination  (sigRA  <  1000  and  sigDec  <  1000); 
see  Munn  et  al.  (2004)  for  details. 


4.  RESULTS 

The  SLoMaSS  observations  were  used  to  investigate  kine¬ 
matics,  colors,  chromospheric  activity,  and  metallicity  in  the  thin 
and  thick  disk  populations. 

We  used  two  methods  in  our  analysis.  First,  we  examined  the 
velocity  distributions  with  no  assumptions  regarding  the  parent 
population  of  a  given  star  (see§4.1).  These  kinematic  results  were 
then  used  to  test  a  leading  Galactic  model  (§  4.2).  We  also  used 
the  space  motion  of  each  star  to  assign  it  to  the  thin  or  thick  disk 
population,  and  investigated  the  color,  metallicity,  and  activity 
differences  between  the  two  populations  (§  4.3). 

4.1.  Kinematics 

The  sample  was  binned  in  100  pc  increments  of  vertical  dis¬ 
tance  from  the  Galactic  plane  and  probability  plots  of  UVW  ve¬ 
locities  were  constructed  (Lutz  &  Upgren  1980;  Reid  et  al.  1995, 
2002).  These  diagrams  plot  the  cumulative  probability  distribu¬ 
tion  in  units  of  the  standard  deviation  of  the  distribution.  Hence, 
a  single  Gaussian  distribution  will  appear  as  a  straight  line  with 
a  slope  corresponding  to  the  standard  deviation  of  the  sample 
and  a  y-intercept  equal  to  the  median  of  the  distribution.  Non- 
Gaussian  distributions  will  significantly  deviate  from  a  straight 
line.  An  example  of  velocity  distributions  and  their  probability 
plots  are  shown  in  Figure  3.  An  advantage  to  this  method  is  its 
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Fig.  5. — Velocity  dispersion  as  a  function  of  vertical  distance  from  the  plane  in  100  pc  bins.  The  open  squares  represent  the  kinematically  colder,  core  component, 
and  the  crosses  represent  the  high-dispersion,  wing  component.  Errors  are  derived  from  the  linear  fits  to  the  probability  plots.  Note  the  smooth  trend  in  the  low- 
dispersion  component  as  a  function  of  height  from  the  plane. 


immunity  to  outliers  and  to  binning  effects  from  poorly  popu¬ 
lated  histograms. 

For  each  distance  bin  in  SLoMaSS,  the  probability  plots  are 
well  fit  with  two  lines:  a  low-dispersion,  kinematically  colder 
“core”  component  (<|1  cr|)  and  a  high-dispersion  “wing”  com¬ 
ponent  (>|1  tr|).  Shown  in  Figure  4  is  an  example  of  our  anal¬ 
ysis,  illustrating  linear  fits  to  the  core  and  wing  distributions.  At 
larger  Galactic  heights,  the  wing  component  traces  the  in  situ 
thick  disk  population. 

The  resulting  dispersions  are  shown  as  a  function  of  distance 
from  the  Galactic  plane  in  Figure  5.  The  low-dispersion  core 
component  is  well  behaved,  smoothly  increasing  with  increasing 
height.  The  high-dispersion  wing  component  is  subject  to  larger 
scatter,  but  generally  increases  with  height  from  the  Galactic  plane. 
Our  results  are  summarized  in  Table  1 .  In  order  to  facilitate  com¬ 
parison  to  the  previous  results  included  in  the  table,  we  report  the 
dispersions  for  several  distance  bins,  as  well  as  the  entire  sample. 

4.2.  Galactic  Dynamical  Models 
4.2.1.  Besangon  Model:  Introduction 

Our  sample  is  well  suited  to  rigorously  test  current  Galactic 
models.  We  chose  the  Besanpon  model  (Robin  et  al.  2003)  as  a 
fiducial,  since  it  simulates  both  star  counts  and  kinematics.  This 
model  is  constructed  on  several  assumptions  and  empirical  con¬ 


straints,  with  the  goal  of  reproducing  the  stellar  content  of  the 
Milky  Way.  Four  populations  comprise  the  model:  the  thin  and 
thick  disks,  bulge,  and  halo.  For  each  population,  a  star  fonnation 
history,  initial  mass  function,  density  law,  and  age  are  imposed.  In 
addition,  an  age-metallicity  relation  is  employed  for  each  popu¬ 
lation,  with  a  Gaussian  dispersion  about  the  mean  metallicity  of 
each  component  (see  Tables  1  and  3  of  Robin  et  al.  2003). 

Kinematically,  the  thin  disk  is  composed  of  seven  groups  of 
different  ages  (from  0.0  to  10  Gyr),  each  being  isothermal  except 
for  the  youngest  (0.0-0.15  Gyr).  The  thick  disk  is  composed  of 
an  1 1  Gyr  old  population,  with  a  velocity  ellipsoid  based  on  the 
measurements  of  Ojha  et  al.  (1996,  1999).  An  age-velocity  dis¬ 
persion  relation  from  Gomez  et  al.  (1997)  is  imposed  on  each 
component,  and  the  model  is  then  allowed  to  self-consistently 
evolve  to  the  present  day.  This  self-consistency  is  achieved  with 
the  method  described  in  Bienayme  et  al.  (1987).  Stellar  popula¬ 
tions  are  formed  according  to  their  appropriate  density  laws  and 
introduced  with  an  initial  velocity  dispersion.  The  mass  density 
of  the  stars  is  summed  in  a  column  of  unit  volume  centered  on  the 
Sun  along  with  dark  matter  halo  and  interstellar  material  contri¬ 
butions,  and  the  potential  is  computed  using  the  Poisson  equa¬ 
tion.  Stars  are  then  evolved  using  the  collisionless  Boltzmann 
equation  in  this  new  potential,  and  redistributed  in  the  z-direction, 
as  complete  orbital  evolution  is  not  included  in  this  model.  The 
process  is  iterated  until  the  potential  and  scale  heights  of  the  disk 
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Fig.  6. — Shown  are  the  R  —  I  and  distance  distributions  of  the  stars  in  SLoMaSS 
( solid  line )  and  a  sampled  model  (dashed  line).  Note  the  agreement  between  the 
two  data  sets,  indicating  that  we  are  sampling  each  instance  of  the  Besarujon  model 
in  a  manner  consistent  with  the  spectroscopic  observations  in  SLoMaSS. 


populations  converge  at  the  1%  level.  The  local-mass  density, 
which  was  empirically  determined  by  Creze  et  al.  (1998),  also 
imposes  a  constraint  on  the  initial  Galactic  potential.  We  note 
that  imposing  an  age-velocity  dispersion  relation  on  the  model 
dominates  the  kinematics,  and  that  there  are  significant  uncer¬ 
tainties  in  these  relations.  For  example,  the  age-velocity  disper¬ 
sions  used  by  Rocha-Pinto  et  al.  (2004  )  differ  by  factors  of  ~2  for 
the  oldest  stars,  compared  to  the  relations  in  Gomez  et  al.  (1997). 
These  discrepancies  are  primarily  derived  from  the  difficulty  in 
determining  ages  of  field  stars,  with  systematic  differences  imposed 
by  various  methods  (i.e.,  isochrone  fitting  vs.  chromospheric  ages). 

Intrinsic  properties  such  as  age,  mass,  luminosity,  metallicity, 
position,  and  velocity  are  available  for  each  star  in  the  simula¬ 
tion.  Observable  properties,  such  as  colors,  proper  motions,  and 
radial  velocities  are  also  reported.  To  determine  colors,  the  model 
uses  the  Lejeune  et  al.  (1998)  database,  which  employs  adjusted 
stellar  synthetic  spectra  that  attempt  to  match  empirical  color- 
temperature  relations. 

4.2.2.  Comparison  to  SDSS 

We  queried  the  Besancon  Web  site,10  generating  a  suite  of 
synthetic  data  sets  along  the  appropriate  Galactic  sight  line  of 


10  See  http://physique.obs-besancon.fr/modele/. 
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Fig.  7. — Probability  plots  of  the  U,  V ,  and  W  velocities  lor  pure  thin  disk  {shaded  line)  and  thin  +  thick  models  (solid  line)  in  the  1200  pc  <  z  <  1300  pc  bin.  This  height 
bin  corresponds  to  the  SLoMaSS  data  shown  in  Fig.  4.  Note  that  when  the  thick  disk  is  added  to  the  distribution,  the  overall  slope  increases,  and  the  wing  components  are 
enhanced,  as  seen  in  the  W probability  plot.  The  combined  thin  and  thick  disk  velocity  distributions  are  also  necessary  to  match  the  two-component  structure  seen  in  the 
SLoMaSS  probability  plots. 
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Fig.  8. — Thin  disk  (top)  and  thick  disk  ( bottom )  velocity  dispersions  for  SLoMaSS  data  ( open  squares)  and  Besanqon  model  predictions  (crosses).  The  model  re¬ 
produces  general  trends,  such  as  increased  dispersion  with  distance  from  the  Galactic  plane.  The  thin  disk  W  velocity  dispersions  are  systematically  low,  as  are  the  thick 
disk  V dispersions  at  large  heights.  This  is  probably  due  to  the  assumed  age-velocity  dispersion  relation,  as  explained  in  §  4.2.4. 


SLoMaSS  (/  ~  105°,  b  ~  —62°)  with  proper  magnitude  limits 
and  error  characteristics.  We  included  both  the  thin  and  thick 
disks  in  the  model  inputs.  The  bulge  and  halo  components  were 
excluded,  since  SLoMaSS  points  away  from  the  bulge,  and  with 
a  maximum  distance  of  ~2000  pc  we  expect  halo  contamination 
to  be  minimal.  A  total  of  25  models  were  generated  from  iden¬ 
tical  input,  in  order  to  minimize  Poisson  noise. 

As  a  consistency  test  of  the  Besanqon  model,  we  compared 
the  model  star  counts  to  those  obtained  with  SDSS  survey  pho¬ 
tometry  in  the  SLoMaSS  field.  We  queried  the  SDSS  Catalog 
Archive  Server11  for  good  point-spread  photometry12  of  stars 
along  the  appropriate  sight  line.  This  query  resulted  in  46,730 
stellar  targets.  In  the  models,  the  mean  star  count  was  46,991 
with  a  standard  deviation  of  137.  This  excellent  agreement 
should  be  viewed  cautiously,  as  only  rough  magnitude  cuts  were 
imposed  on  the  Besan$on  models,  and  we  did  not  attempt  to 
model  observational  problems  that  would  affect  SDSS  star  counts 
(Vanden  Berk  et  al.  2005),  such  as  cosmic  rays  or  diffraction 
spikes  near  bright  stars.  However,  the  test  inspires  confidence 


11  See  http://cas.sdss.org/dr5/en/. 

12  Specifically,  we  required  the  following  flags:  a  detection  in  BINNED  1  and  no 
EDGE,  NOPROFILE,  PEAKCENTER,  NOTCHECKED,  PSF_FLUX_INTERP, 
SATURATED,  BAD_COUNTS_ERROR,  DEBLEND- NOPEAK,  CHILD, 
BLENDED,  INTERP-  CENTER,  or  COSMIC-RAY  flags  set  (Stoughton  et  al. 
2002). 


that  the  Besanqon  model  adequately  represents  the  observed  SDSS 
star  counts. 

4.2.3.  Comparison  to  SLoMaSS 

Each  of  the  25  models  was  subsampled  to  reproduce  the 
color  and  distance  distributions  observed  in  SLoMaSS.  The  two- 
dimensional  color-distance  density  distribution  was  computed 
for  SLoMaSS,  and  color-distance  pairs  were  randomly  drawn 
according  to  this  normalized  density  distribution.  If  the  color- 
distance  pair  was  found  in  the  Besanijon  model,  then  it  was  kept. 
This  Monte  Carlo  sampling  continued  until  there  were  7398 
color-distance  pairs,  matching  the  number  of  stars  in  SLoMaSS. 
This  sampling  forced  each  model  to  simulate  the  properties  of 
the  SLoMaSS  observations  when  compared  to  the  complete 
SDSS  photometry.  Since  SDSS  r  —  i  colors  are  not  available  for 
the  Besan9on  model,  the  SLoMaSS  colors  were  transformed  to 
R  —  I  using  the  relations  of  Davenport  et  al.  (2006).  Shown  in 
Figure  6  are  the  R  —  I  and  distance  distributions  for  one  instance 
of  the  model,  along  with  those  from  SLoMaSS.  Thus,  each  model 
is  “observed”  in  a  manner  consistent  with  the  spectroscopic  ob¬ 
servations  of  the  SDSS  field  photometry. 

4.2.4.  Kinematic  Comparisons 

Using  the  subsampled  model  data,  we  compared  the  kinematic 
predictions  of  the  Besanqon  models  to  the  UVW  velocities  of 
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Fig.  9. — Thin  disk  (top)  and  thick  disk  (bottom)  mean  velocities  for  SLoMaSS  data  (open  squares)  and  Besangon  model  predictions  (crosses).  The  thin  disk  mean  V 
velocities  predicted  by  the  Besangon  models  are  systematically  low,  while  the  observed  thick  disk  mean  V  velocities  diverge  at  large  Galactic  heights. 


SLoMaSS.  Using  the  method  described  in  §  4.1,  we  constructed 
two  sets  of  probability  plots  for  each  model:  one  composed  of  the 
thin  and  thick  disk  stars  and  one  measured  solely  for  the  thin 
disk.  The  reason  for  this  separation  was  twofold.  The  first  was  to 
isolate  any  systematics  between  the  models  thin  and  thick  disk 
predictions.  In  addition,  this  allowed  for  direct  testing  comparison 
of  the  thin  disk  predictions,  since  most  of  these  stars  lie  on  the 
inner  core  region  of  the  probability  plot.  That  is,  if  the  velocity 
predictions  are  correct  for  the  thin  disk,  the  overall  slope  of  the 
isolated  thin  disk  models  should  roughly  match  the  core  slope 
measured  from  SLoMaSS. 

In  Figure  7,  an  illustrative  example  of  this  analysis  is  shown. 
The  main  effect  of  adding  the  thick  disk  component  is  to  increase 
the  slope  of  both  the  core  and  wing  components.  In  addition,  the 
wing  component  is  enhanced,  as  seen  in  the  W  velocity  distribu¬ 
tion,  which  is  expected  from  the  addition  of  a  high-dispersion 
population.  It  is  clear  from  comparing  to  Figure  4  that  the  com¬ 
bination  of  the  thin  and  thick  disk  model  predictions  are  neces¬ 
sary  to  simulate  the  structure  seen  in  the  SLoMaSS  data. 

Following  the  method  explained  in  §  4.1  we  measured  the 
slopes  of  the  core  and  wing  components  of  the  each  model  as 
a  function  of  distance  from  the  Galactic  plane.  The  results  of 
this  analysis,  compared  to  the  SLoMaSS  results  are  shown  in 
Figures  8  and  9.  The  SLoMaSS  results  ( open  squares),  which 
are  shown  in  Figure  5,  are  compared  to  the  average  Besangon 


prediction  ( crosses )  for  the  thin  (top)  and  thick  disk  (bottom). 
While  the  model  does  well  in  predicting  general  trends,  there 
are  clearly  some  systematics.  The  model  predictions  for  the  o> 
thin  disk  velocity  dispersions  are  systematically  low,  suggest¬ 
ing  there  may  be  flaws  in  the  method  used  to  compute  these 
motions.  In  addition,  the  model  underestimates  the  thick  disk 
cry  dispersions  at  large  Galactic  heights.  As  described  above 
(§  4.2),  the  assumed  age-velocity  dispersions  relations,  which 
dominate  the  predicted  kinematic  structure,  are  uncertain  and 
may  contribute  to  this  disagreement.  This  speculation  is  sup¬ 
ported  by  the  results  summarized  in  Table  1,  which  demon¬ 
strate  that  most  previous  surveys  (as  well  as  our  own)  have 
measured  velocity  dispersions  significantly  higher  than  those 
predicted  from  the  model.  The  mean  velocity  distributions  for 
both  SLoMaSS  and  the  Besanc^on  model  are  shown  in  Figure  9. 
Again,  the  model  in  general  performs  well,  but  there  are  evi¬ 
dent  systematic  differences.  The  mean  V  velocities  predicted 
for  the  thin  disk  by  the  model  are  systematically  high,  and  the 
thick  disk  V velocities  diverge  at  large  Galactic  heights.  The  first 
difference  may  be  attributed  to  our  chosen  solar  motion  values 
(10,  5,  and  7  km  s_1;  Dehnen  &  Binney  1998).  A  larger  adopted 
value  of  Vq,  such  as  the  classic  value  of  12  km  s_1  (Delhaye 
1965),  would  move  the  mean  velocities  toward  agreement.  The 
discrepancy  in  the  thick  disk  V velocities  is  probably  due  to  small 
number  statistics,  as  seen  in  Figure  2. 
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Fig.  10. — Metallicity-sensitive  ratio  of  (CaH2  +  CaH3)/Ti05  vs.  spectral  type 
for  the  thin  (open  squares)  and  thick  ( crosses )  disk  populations.  Higher  ratio  val¬ 
ues  indicate  a  higher  metallicity  (Burgasser  &  Kirkpatrick  2006).  Within  the  error 
bars,  the  thin  and  thick  disk  populations  exhibit  very  similar  behavior,  indicating 
that  the  sample  may  not  be  probing  a  large  spread  (  ^  1  dex)  in  metallicity. 


Spectral  Type 

Fig.  1 1 . — The  ugriz  colors  as  a  function  of  spectral  type  for  the  thin  (open 
squares)  and  thick  ( crosses )  disk  populations.  The  colors  of  both  populations 
agree  within  the  errors,  suggesting  that  the  M  dwarfs  in  SLoMaSS  do  not  cover  a 
large  metallicity  range. 


4.3.  Differences  between  the  Thin  and  Thick  Disks 

In  order  to  examine  observable  differences  between  M  dwarf 
members  of  the  thin  and  thick  disks,  we  kinematically  sepa¬ 
rated  the  sample  using  the  method  of  Bensby  et  al.  (2003).  This 
technique  selects  outliers  in  the  wings  of  the  three-dimensional 
Gaussian  velocity  distribution,  and  computes  the  probability  of 
these  stars  belonging  to  the  thick  disk.  The  thin  and  thick  disk 
populations  were  characterized  by  the  velocity  dispersions  in 
Bensby  et  al.  (2003).  In  order  to  minimize  systematics,  the  UVW 
velocities  were  recomputed  using  distances  detennined  with  the 
spectroscopic  parallax  relations  of  Hawley  et  al.  (2002)  as  de¬ 
scribed  in  §  3.3.  Thus,  UVW  should  not  a  priori  vary  systemati¬ 
cally  with  color.  This  also  limits  the  analysis  to  the  6577  M  dwarfs 
in  the  sample. 

4.3.1.  Metallicity 

The  exact  formation  mechanism  of  the  thick  disk  is  uncertain 
(see  Majewski  1993  and  references  therein),  but  there  is  evi¬ 
dence  that  the  mean  metallicity  of  the  thick  disk  is  lower  than  that 
of  the  thin  disk  (Reid  &  Majewski  1993;  Chiba  &  Beers  2000; 
Bensby  et  al.  2003).  In  addition,  differences  in  the  a -element 
distributions  have  been  shown  to  be  distinct  (Fuhrmann  1998; 
Feltzing  et  al.  2003;  Bensby  et  al.  2003).  While  direct  metallicity 
determinations  of  M  dwarfs  are  difficult  (Woolf  &  Wallerstein 
2006),  proxies  of  metallicity  have  been  employed  (Gizis  1997; 
Lepine  et  al.  2003;  Burgasser  &  Kirkpatrick  2006).  These  previ¬ 
ous  studies  have  used  the  low-resolution  molecular  band  in¬ 
dices  (CaHl,  CaH2,  CaH3,  and  Ti05)  defined  in  Reid  et  al. 
( 1995)  to  roughly  discriminate  between  solar-metallicity,  sub¬ 
dwarf  ([m/H]  ~  —1.2),  and  extreme-subdwarf  ([m/H]  ~  —2) 
populations.  Adapting  the  methods  of  Lepine  et  al.  (2003)  and 
Burgasser  &  Kirkpatrick  (2006)  we  computed  the  ratio  (CaH2  + 
CaH3)/Ti05  for  each  star  in  the  sample,  which  varies  such  that 


a  larger  value  indicates  a  higher  metallicity.  The  mean  ratio  for 
each  spectral  type  for  the  thin  and  thick  disk  populations  is  shown 
in  Figure  10.  The  two  populations  do  not  separate  within  the 
errors,  suggesting  that  the  observed  metallicity  distributions  do 
not  differ  greatly. 

4.3.2.  Color  Gradients 

A  lower  mean  metallicity  could  also  manifest  itself  as  a  color 
shift  at  a  given  spectral  type.  Specifically,  West  et  al.  (2004) 
showed  that  at  a  given  r  —  i  or  i  —  z  color,  low-metallicity  sub¬ 
dwarfs  ([m/H]  ~  — 1.2)  were  redder  in  g  —  r.  After  SLoMaSS 
was  kinematically  separated,  the  mean  u  —  g,g  —  r,  r  —  i,  and 
i  —  z  colors  of  the  thin  and  thick  disk  stars  were  computed  at 
each  spectral  subtype,  as  shown  in  Figure  11.  There  are  no  sig¬ 
nificant  color  differences  between  the  thin  and  thick  disks  at  a 
given  spectral  type.  This  further  suggests  that  the  SLoMaSS  stars 
are  not  probing  a  large  spread  in  metallicity. 

4.3.3.  Chromospheric  Activity 

Finally,  if  the  thick  disk  is  an  older  system  (as  suggested  by 
its  hotter  kinematics),  it  should  possess  a  lower  fraction  of  active 
stars  (West  et  al.  2004,  2006).  The  chromospheric  activity  time- 
scale  varies  with  mass,  such  that  higher  mass  stars  lose  their  ac¬ 
tivity  sooner  (after  ~  1  Gyr)  than  their  low-mass  counterparts 
('-MO  Gyr).  To  observationally  quantify  chromospheric  activity, 
the  Ha  equivalent  width  (EW)  is  measured  for  each  M  dwarf. 
We  employed  the  technique  introduced  in  West  et  al.  (2004)  and 
described  in  Bochanski  et  al.  (2007),  selecting  chromospherically 
active  and  inactive  stars  at  each  spectral  type.  To  be  classified  as 
active,  a  star  must  have  an  EW  of  1  A,  and  pass  additional  signal- 
to-noise  ratio  and  error  tests  described  in  West  et  al.  (2004). 
Figure  12  shows  the  active  fraction  of  stars  as  a  function  of 
spectral  type  for  both  disk  populations.  While  most  early-type 
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Fig.  12. — Active  fraction  of  stars  as  a  function  of  spectral  type  for  the  thin 
(open  squares )  and  thick  (crosses)  disk  populations.  An  older  population  would 
show  a  lower  active  fraction,  as  suggested  by  the  thick  disk  stars  in  the  M5  bin. 
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further  investigation  into  modeling  the  kinematics  of  the  thin  and 
thick  disks  is  necessary.  In  particular,  the  age-velocity  dispersion 
plays  an  important  role  in  the  kinematic  predictions,  and  a  more 
exact  definition  of  this  relation  is  required. 

SLoMaSS  was  divided  kinematically,  assigning  membership 
of  each  star  to  either  the  thin  or  thick  disk.  We  inspected  the  two 
populations  as  functions  of  spectral  type,  searching  for  observ¬ 
able  differences  in  the  metallicity,  colors  and  chromospheric 
activity.  While  there  was  little  observed  difference  between  the 
metallicity  and  color  distributions,  the  activity  fraction  distribu¬ 
tion  suggests  that  the  thick  disk  displays  an  activity  level  con¬ 
sistent  with  an  older  population,  although  this  result  is  purely 
suggestive  and  needs  to  be  re-examined  with  later  spectral  types. 
The  lack  of  a  strong  observable  signal  may  have  several  causes. 
Primarily,  kinematic  separation  of  populations  is  not  perfect,  as 
stars  with  extreme  kinematics  in  one  population  (i.e.,  the  thin 
disk)  can  masquerade  as  members  of  the  other  population.  This 
would  dilute  observable  differences,  and  in  our  case,  the  more 
numerous  thin  disk  population  may  be  polluting  the  thick  disk 
sample,  despite  our  best  efforts  to  minimize  this  effect.  Second, 
the  intrinsic  properties  (i.e.,  metallicity)  of  the  stars  that  comprise 
the  thin  and  thick  disks  are  drawn  from  overlapping  distribu¬ 
tions.  Stars  that  are  at  the  edges  of  these  distributions  would  also 
blur  the  distinction  among  observable  properties,  clouding  the 
best  efforts  of  the  observer.  Further  investigations  using  stellar 
tracers  along  the  entire  main  sequence  should  alleviate  some  of 
these  issues. 


M  dwarfs  (M0-M3)  lose  their  activity  quite  rapidly,  the  later 
types  (M5)  possess  smaller  active  fractions  in  the  thick  disk,  and 
exhibit  the  behavior  expected  of  older  systems.  However,  these 
results  are  only  suggestive,  as  populations  older  than  ~4  Gyr 
require  M  dwarfs  with  types  later  than  M5  in  order  to  be  effective 
probes  of  the  age  (A.  A.  West  et  al.  2008,  in  preparation). 

5.  CONCLUSIONS 

The  Milky  Way  (along  with  other  spiral  galaxies)  is  evidently 
a  composite  of  a  few  major  smooth  components  (the  thin  and 
thick  disks  and  the  halo)  and  many  smaller  structures,  such  as  the 
tidal  debris  streams  from  the  Sagittarius  dwarf  galaxy.  Using 
spectra,  proper  motions,  and  photometry  along  one  Galactic 
sight  line,  we  studied  the  kinematic  and  structural  distributions 
of  the  smooth  thin  and  thick  disks  as  a  function  of  vertical  distance 
from  the  plane.  We  fit  two-component  Gaussians  to  the  UVW 
distributions  as  a  function  of  height,  placing  new  constraints  on 
the  thin  and  thick  disk  velocity  dispersions. 

This  sample  was  then  employed  to  test  the  predictions  of  cur¬ 
rent  Galactic  models.  The  Besantjon  model  was  chosen  for  com¬ 
parison,  since  it  is  widely  accepted  as  a  standard  and  models  the 
kinematics  of  stellar  populations.  We  generated  a  suite  of  25  mod¬ 
els,  and  each  model  was  sampled  to  be  consistent  with  the  col¬ 
ors  and  distance  distribution  of  SLoMaSS.  The  bulk  kinematic 
properties  of  the  data  and  model  were  compared.  They  agree  to 
~  1 0  km  s  1 ,  placing  a  limit  on  the  accuracy  of  model  predic¬ 
tions.  However,  aw  is  poorly  fit  by  the  model,  suggesting  that 
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